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Constraints on big bang nucleosynthesis (BBN) and on cosmological pa- 
rameters from conflicting deuterium observations in different high red-shift 
QSO systems are discussed. The high deuterium observations by Carswell 
et al, Songaila et at, and Rugers Sz Hogan is consistent with 4 He and 7 Li 
observations and Standard BBN (N u =3) and allows N v < 3.6 at 95% C.L., 
but is inconsistent with local observations of D and He in the context of 
conventional theories of stellar and Galactic evolution. In contrast, the low 
^ ■ deuterium observations by Tytler, Fan & Buries and Buries & Tytler are 

consistent with the constraints from local Galactic observations, but require 
N v = 1.9 ± 0.3 at 68% C.L., excluding Standard BBN at 99.9% C.L., unless 
the systematic uncertainties in the 4 He observations have been underesti- 
mated by a large amount. The high and low primordial deuterium abun- 
dances imply, respectively, Q B h 2 = 0.005 - 0.01 and Q B h 2 = 0.02 - 0.03 
at 95% C.L. When combined with the high baryon fraction inferred from x- 
ray observations of rich clusters, the corresponding total mass densities (for 
50 < H < 90) are ft M = 0.05 - 0.20 and ft M = 0.2 - 0.7, respectively (95% 
C.L.) The range of corresponding to high D is in conflict with dynam- 
ical constraints (f^M > 0.2 — 0.3) and with the shape parameter constraint 
(r = £ly[h = 0.25 ± 0.05) from large scale structure formation in CDM and 
ACDM models. 
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I. INTRODUCTION 



Among the light nuclides synthesized during the early evolution of the universe, deu- 
terium is unique in its sensitivity to the universal density of baryons and in the simplicity of 
its galactic evolution. As gas is incorporated into stars and the heavy elements are synthe- 
sized, D is only destroyed |TJ so that any D abundance inferred from observations provides 
a lower bound to its primordial value. [] Unfortunately, an upper bound to the primordial 
D abundance is more uncertain, depending on the evolutionary history of the matter being 
observed. Thus, although estimates of the D abundance in the presolar nebula and in 
the local interstellar medium (ISM) provide interesting lower bounds to primordial D 
[X 2P > X 2Q = (3.6 ±1.3) x 1(T 5 , X 2P > X 2lSM = (2.2 ± 0.3) x 1(T 5 , where the H mass 
fraction has been taken to be X Q = Xi$ M = 0.70 ± 0.01], upper bounds are more model 
dependent (see, for example, Ref. For this reason, observations of D in (nearly) 

unevolved systems (high red-shift, low metallicity QSO absorbers) have been eagerly antic- 
ipated. If, indeed, A 2 p ~ A 2 qso, then because of the sensitivity of the D abundance to 
the nucleon abundance (rj = n^jn^] the ratio of the present baryon density to the critical 
density is Qp,h 2 = 0.0037?7io, where the Hubble parameter is H = 100/i km/s/Mpc and 
r]i = 10 w r]), a measurement of (D/H)q SO to ~30% accuracy will lead to a determination 
of r) to ~ 20% accuracy. Armed with 77, reasonably accurate predictions of the primordial 
abundances of 3 He, 4 He, and 7 Li will follow (see, for example, Ref. For example, for 
1.5 < i] W < 10, a 20% uncertainty in 7710 will lead to an uncertainty in the predicted 4 He 
mass fraction which ranges from ~ 0.003 (at low 7710) to ~ 0.002 (at high rjio) . Deuterium 
is the ideal baryometer. 

In the last two years, observations of D in high red-shift, low metallicity QSO absorbers 



1 X 2 p > ^20BS; where the D mass fraction is X 2 = 2X nz)/nn; X is the hydrogen mass fraction, 
and n x is the number density for nuclide x; the subscript P is for the primordial abundance. As 
an estimate of the primordial value Xp = 1 — Yp, we will adopt Xp = 0.76 ± 0.01. In this paper 
we quote la uncertainties unless otherwise indicated. 
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have begun to appear in the published literature JT0|-|IJ]. The first observations of D in 
absorption against Q0014+813 [p~0| p~2[] suggested a surprisingly high abundance [for our 
quantitative comparisons we will adopt the recent reanalysis by Rugers and Hogan: D/H 
= (1.9 ± 0.4) x 10~ 4 , X 2 = (2.9 ± 0.6) x 10~ 4 ], roughly an order of magnitude larger 
than the presolar or ISM values (A 2 q S o/A 2ISM ~ 8 ± 3, X 2 q S o/X 20 ~ 13 ± 3). As such 
efficient D destruction in the Galaxy is not expected [PH8|J3|JT5[| it has been suggested that 
the feature identified as D in Q0014+813 might be a hydrogen interloper |16|]. However, 



the Rugers-Hogan reanalysis argues against this possibility. Further, recent papers |17|jl8 
present evidence for D absorption in front of two other QSOs (Q0420-388 and BR1202- 
0725, respectively) which, if the identifications are correct, suggest D/H > 2 x 10~ 5 and 
D/H < 1.5 x 10~ 4 , respectively. Although puzzling from the point of view of chemical 
evolution in the Galaxy, the high D abundance points towards a low baryon density (r]i ~ 2) 
which is consistent with the predicted and observed (inferred) primordial abundances of 4 He 
and 7 Li [|16|,[n|. As we shall see, however, this low baryon density (Q^h 2 ~ 0.007) is in 
conflict with determinations of the total mass density and the baryon fraction inferred from 
x-ray observations of rich clusters. 

In contrast, from recent observations, Tytler, Fan, and Buries [|I3] and Buries and 
Tytler @ derive a low D abundance: (D/H) = [2.3 ± 0.3 (stat) ±0.3 (sys) ] x 10~ 5 



towards the QSO1937-1009 @ and (D/H) = [2.5±g;t (stat) t° i (sys)] xl0~ 5 towards 
the QSO1009+2956 @. We have combined their two results to obtain (D/H) QSO = 
(2.4 ± 0.5) x 10" 5 ; X 2QS o = (3.6 ± 0.8) x 10~ 5 . Although marginally larger than ISM 
deuterium (X 2 q SO /X 2ISM = 1-6 ± 0.4), the low D abundance is not very different from the 
presolar value (X 2 qso/A 20 = 1.0 ± 0.4), suggesting that even though the absorbers are at 
high redshift (z a b s = 3.572 and 2.504) and have very low metallicity (~ 10~ 3 solar), some 
D may have already been destroyed (A 2 p > A 2 q SO ). If indeed A 2P ~ A 2 q SO (no significant 



D destruction), then the problems for BBN identified by Hata et al. ||19|| , which were based 
on A 2 p inferred from solar system observations of D and 3 He, persist. The higher baryon 
density suggested by the low D result is, however, in agreement with the x-ray cluster data 



(but still supports a low density universe). 

It is hoped that future observations of D in other high red-shift, low metallicity QSO 
absorbers will resolve the current dichotomy between the high D result for Q0014+813 
on the one hand 0-0 and the low D results for Q1937-1009 and Q1009+2956 on the 
other [|n],[n]]. Here, we explore the implications for cosmology (the baryon density), for the 
primordial abundances of the other light nuclides ( 4 He and 7 Li), and for particle physics 
(bounds to the effective number of equivalent light neutrinos, N v ), of the high D abundance 
and contrast them with those for the low D abundance. 

II. PRIMORDIAL D AND BBN 

If (D/H)p is fixed, standard big bang nucleosynthesis (SBBN: homogeneous, N v = 3, 
the neutron life time r n = 887 ± 2 s, etc) can be used to predict the primordial abundances 
of the other light nuclides and determine the present baryon density. In Fig. [l] the SBBN 
predicted abundances of 4 He (mass fraction Y), D (y 2 = D/H), and 7 Li (y 7 = 7 Li/H) are 
shown as a function of the nucleon to photon ratio r\ for 1 < 7710 < 10. Convolving the SBBN 
predictions (including uncertainties estimated by the Monte Carlo method of Ref f2"fj|| ) with 
the high D and low D results constrains 77 and leads to predictions of Y and y-j as may 
be seen in Fig. D], Also shown in Fig. |T] are the 68 and 95% C.L. contours for the overlap 
between the inferred primordial abundances of 4 He [Yp = 0.232 ± 0.003 (stat) ±0.005 (sys)] 
and 7 Li [log y-j = —9.8 ± 0.2 (sys) ±0.3 (depletion/creation)] Q and the BBN predictions. 
When using the low D value we must ensure that consistency with the ISM (and solar 
system) value is maintained (X 2P > X 2 i SM ; X 2P > X 2Q ). In Fig. || we show the SBBN 
likelihood distribution (solid curve) for the low D result [\ogy® so = —4.62 ± 0.05 ± 0.06; 
^qso _ ^ 4 -j- 0.5) x 10~ 5 ]. Requiring that the QSO D abundance be no smaller than 



2 The statistical uncertainty in 4 He is assumed to be Gaussian, while the systematic uncertainties 
in 4 He and Li and the uncertainty in 7 Li depletion/creation are treated as flat (top hat) distribu- 
tions. The statistical uncertainty in y-j is small compared to the systematic and depletion/creation 
uncertainties. 
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the ISM D abundance [y^ so > yf M = (1.6 ± 0.2) x 1(T 5 ], modifies the distribution to 
the dotted curve in Fig. 0, very slightly truncating the lower end of the y 2 distribution 
^qso/ism _ ^4 _|_ q ^ x io -5 ]. It is this latter distribution which we will use in our 
comparisons. If, further, we also require that the QSO D abundance exceeds that inferred 
for the presolar nebula [y 2 Q = (2.6 ± 0.9) x 10~ 5 ] [0, the distribution (dashed curve in 
Fig. g) is slightly shifted to higher values [^so/ism/© = (2.6 ± 0.5) x 10~ 5 ]. Given that the 
solar system data may be subject to different systematic errors than those associated with 
the QSO and ISM absorption observations, we will limit our analysis to those which follow 
from the marginally less restrictive QSO/ISM constraint. The resulting SBBN constraints 
on 7), VL^li 2 , N„, Y, and y-j which follow from the high and low-QSO/ISM D abundances are 
summarized in Table | along with, for comparison, the previous Hata et al. results @ which 
utilized solar system D and 3 He abundances. Fig. ^ shows the likelihood distributions for rj 
for high and low primordial deuterium. 

A glance at Fig. [l| reveals the well-known result [21| that high-D is consistent with the 
primordial 4 He abundance inferred from HII region data. In Fig. [| is shown the SBBN 
predicted 4 He mass fractions corresponding to the two deuterium values. It is clear from 
Figs. [I] and |] that, in the absence of large systematic errors in 4 He, the observed 4 He 
abundance favors high-D and is inconsistent with low-D. 

In Fig. [| we compare the SBBN predictions of 7 Li for high and low deuterium with that 
inferred from observations of the Pop II halo stars |22|-^5|. As is clear from Figs. [I] and |5|, 
and Table [IJ consistency with lithium is achieved for both high and low D. Notice, however, 
that while the high-D (low rj) overlap with SBBN bounds stellar destruction/dilution of 7 Li 
(y7 < 3.0 x 10~ 10 ; log y-j < —9.5), the low-D (high rj) overlap actually requires some modest 
destruction/dilution (3.0 < 10 10 ?/ 7 < 7.8; -9.5 < log y 7 < -9.1). 

Finally, we turn to 3 He, whose post-BBN evolution is model (galactic chemical evolution) 
dependent. However, since production of 3 He by low mass stars can only increase the 3 He 
abundance, observations of D and 3 He constrain the primordial D and 3 He abundances 
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7|j3|j9|] • For any chemical evolution history the observed and primordial abundances 
of D and 3 He may be related through one parameter g 3 , the effective 3 He stellar survival 
fraction, which contains all the stellar and galactic evolution uncertainties [p6|,|3|,|9f . While 



for a single generation of stars g% > 0.25 [28, 29 1 , many specific evolution models suggest 



93 > 0.5 . Following Ref. || we show in Fig. ^| the allowed regions in the Usp-yzp 

plane inferred from SBBN and the high and low deuterium abundances. Although the low-D 
data is entirely consistent with the galactic evolution of D and 3 He, the high-D data requires 
a surprisingly small value of g% (< 0.10 at 95% C.L.) for consistency. Indeed, the high-D 
data suggests that more than 90% of the present ISM has been cycled through stars (since 
D would have to have been destroyed by a factor of ~ 10). With such efficient processing of 
gas through stars, the low metallicity of the ISM is a challenge to galactic chemical evolution 
models HfFl. 



III. SBBN AND N v 



From the discussion above it is clear that high primordial D is entirely consistent with 
the predictions of SBBN and the observed abundances of 4 He and 7 Li. For low primordial D 
there is a significant tension between the predictions of SBBN and the inferred primordial 
abundance of 4 He. If we allow N v , the equivalent number of light neutrinos, to depart 
from the SBBN value N v = 3, we may use the combined D, 4 He, and 7 Li data to find the 
best N u . Fig. [7] shows the N u likelihood distributions for high and low (QSO/ISM) D; for 
comparison we also show the distribution derived from solar system D and 3 He with = 
0.25 — 0.50 |J. Of course, it is always possible that Yp, inferred from nearly primordial (low 
metallicity) extragalactic HII regions |HH 3~3"|1, has been underestimated due to systematic 
errors. Although such uncertainties (ionization corrections, collisional excitation corrections, 
corrections for dust, corrections for stellar absorption, etc.) could either decrease or increase 
the inferred value (Yp = 0.232 ± 0.003), recent work has emphasized those corrections which 

might increase Y P US If we write Yp = 11232 ± °' 003 + AYs y s ( where AY ^ ^ °)> then 
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there is a direct relation between Ay sys and N u , which we show for high and low deuterium 
in Fig. g. For AF sys < 0.009, high-D and SBBN are consistent (N u = 3), but, if AF sys is 
larger, N v > 3 would be required. In contrast, SBBN and low-D are inconsistent unless 



AY sys > 0.011. 



IV. O b AND n 



M 



Either choice of high or low deuterium leads, through SBBN, to reasonably tight con- 
straints on the present ratio of nucleons to photons (see Fig. |I| and Table |), thus bounding 
the present universal density of baryons pe- Comparing pe to the critical density p c , we 
have 

tt B = — = 3.66 x 10- 3 r] W h- 2 . (1) 

Pc 

In Fig. § we show the Qb vs. Hq relation, where the two bands correspond to the 68 and 
95% CL. ranges for r]io allowed by the QSO D abundances (see Table H). Also shown in 
Fig. |9| is an estimate of the luminous baryons identified by observations in the radio, optical, 
ultra-violet, and x-ray parts of the spectrum |j36| , 



Q LVM = 0.004 + 0.0007/i~ 3/2 . (2) 

Over the entire range of H , Qb > ^lum, suggesting the presence of dark baryons. As a 
minimum estimate of the total density (baryons plus non-baryonic dark matter) inferred 
from the dynamics of groups, clusters, etc., we have adopted ^dyn > 0.2 |37|-|39| (although 
others have suggested Qbyn ~ 0.3 [fTO] ). As can be seen from Fig. |, unless H is very small 
(and D is very low), fi DY N > ^b, providing support for the presence of non-baryonic dark 
matter. 

X-ray emission from the hot (baryonic) gas in rich clusters of galaxies offers a valuable 
probe of the fraction of the total mass in the universe contributed by baryons. Although 
relatively rare, such large mass concentrations are expected to provide a fair sample of 
/ B = ^b/^m, where Qm is the total matter density parameter [^l|-|43|]. For clusters, /b = 



Mb/Mtot > M-hq/Mtot, where Mhg is the mass of the x-ray emitting hot gas in the 
cluster and Mtot is the total mass which determines the cluster binding. It is conventional 
to write /hg = /so^so^ 2 , where h$o = H /50 km/s/Mpc, so that 

ft M /^ 2 < 0.0146W/50- (3) 

The inequality in (3) arises from the neglect of the baryons in the galaxies of the cluster; their 
inclusion would reduce the upper bound on Qm by ^ 5-20% (depending on H ). However, 
since the presence of other (dark) baryons (e.g., Machos) cannot be excluded observationally 
and may be large j44], the cluster data is best utilized to provide an upper bound to Qm- 



Thus, in our subsequent analysis we shall employ x-ray data and BBN to evaluate the right 
hand side of Eqn. [3] which we will use to provide an upper bound to Qm (as a function of 
Ho). 

The surprise provided by x-ray observations of rich clusters has been the relatively large 
baryon fraction (/ 50 > 0.1 — 0.2) which, when coupled to the relatively low upper bound 
on r] w from BBN, has led to the "X-Ray Cluster Baryon Catastrophe" PT|,i2"|,i5|: < 1 
unless H is very small. 



Following the recent analysis of Evrard et al. and Evrard we adopt /j 



/ 50 



0.20±0.03, and use this and the bounds on 7710 from high and low D (see Table [j]) to constrain 



the VLm vs. H relation (Eqn. y) in Fig. [L0| We have allowed H to remain unconstrained 
although we believe that recent data suggest H Q = 70 ± 15 km/s/Mpc (h = 0.7 ± 0.15 and 
h 5Q = 1.4 ± 0.3). The x-ray cluster constraints require low Q M , excluding the preferred 
Einstein-de Sitter value of Qm — 1 unless /50 and/or H is much smaller than data indicate. 

Further evidence for low in the context of cold dark matter models (CDM) comes 
from large scale structure constraints on the shape parameter (see, for example, |P7[): T = 
f2 M /i = 0.25 ± 0.05. Since the popular inflationary paradigm suggests that the 3-space 
curvature may vanish, evidence in favor of low Qm has led to consideration of an alternative 
cosmology with a non-vanishing cosmological constant (A) such that ^tot = + — 1 
(see, for example, ||45|| ). Such A cold dark matter models (ACDM) provide the additional 
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benefit of helping to resolve the "age problem" (for the same value of Qm, and fixed H , a 
ACDM universe is older than the corresponding CDM universe). In Figs. and ^ we show 
(at 68% C.L.) the regions in the Hq-Qm plane consistent with the x-ray/BBN constraints 
(for high and low D) and with the shape parameter (r). Also shown is the Qm~Hq relation 
for two choices of the present age of the universe (to = 12 and 15 Gyr) along with the 
dynamically inferred lower bound to the mass density (IIdyn)- I n both cases (CDM and 
ACDM), the low-D, high-?7io choice is preferred over the high-D, low-?yio result. 

A third popular cosmology is the mixed, hot plus cold dark matter model (HCDM; see for 
example, f48H5"0|). In its standard version it is assumed that Om = 1, but that ~ 20 — 30% 
of Qm is in hot dark matter (e.g., neutrinos with mass of 1 - 10 eV) which is relatively 
unclustered on large scales. The shape parameter constraint is not relevant to constraining 
the HCDM model, but the requirement that ^tot — 1 exacerbates the x-ray cluster baryon 
catastrophe pl|- ^ , ^5|1 and the age problem. Even if all HDM could be excluded from x-ray 



clusters which seems unlikely |52| , ~ 1 — ^hdm- Coupled with the upper bound 
from the x-ray cluster data, Qm < 0.7, this requires Qhdm ~ 0.3, nearly closing the preferred 
window (0.2 < n HDM ^ 0.3) on HCDM models. 

V. DISCUSSION 



A determination of the deuterium abundance in a nearly uncontaminated environment 
such as that provided by high redshift, low metallicity QSO absorption clouds could be a 
key to testing the consistency of primordial nucleosynthesis in the standard, hot, big bang 
cosmology, to pinning down the universal density of baryons, and to constraining physics 
beyond the standard model of particle physics. Such data is beginning to be acquired but, 
at present, the observational situation is in conflict. On the one hand there is evidence in 
favor of high D QlllSBIIllOl : (D/H)~ 2 x 10" 4 . In contrast, Tytler, Fan, & Buries g3 



and Buries & Tytler |14| find evidence for low D: (D/H)~ 2 x 10 5 . If the former, high-D 
values are correct, it is surprising that Tytler, Fan, & Buries and Buries & Tytler fail to 
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find such a large abundance in their high redshift (z = 3.57 and 2.50), very low metallicity 
(Z/Z Q ~ 10 -3 ) absorbers; high z and low Z argue against an order of magnitude destruction 
of primordial D. If, instead, the low D result is correct, such weak D-absorption might 
often go unnoticed and the high-D cases might be accidental interlopers. Based on velocity 
information Rugers and Hogan |T2| argue against this possibility which, if more high-D cases 



are found, will become increasingly unlikely. Presumably, the present confused situation 
will be clarified by the acquisition of more data. Here, we have considered separately the 
consequences for cosmology and particle physics of the high-D and low-D data. 

For the high-D case, SBBN (N u = 3) is consistent with the inferred primordial abun- 
dances of D, 4 He, and 7 Li provided that the baryon density is small (see Table | and Figs. [I], 
|3|, [|, [5|, |7|, and H). However, for consistency with the solar system and/or present interstellar 
D abundances, such a large primordial D abundance requires very efficient D destruction. 
The low baryon density which corresponds to high-D still leaves room for dark baryons and 
reinforces the case for non-baryonic dark matter (see Fig. |9[). However, when combined with 
the x-ray cluster data, low Q B and high / H g suggest a very low density universe (Om < 0.21 
for H Q > 50 km/s/Mpc; see Fig. ITDp . The conflict between the upper bound on Q M and the 



evidence for a lower bound ^dyn ~ 0.2 — 0.3 argue against high-D and low r] 10 (see Figs. [TT] 



and |T|) . Such a low value for is also in conflict with the constraint from the shape 



parameter V (see Figs. |ll] and fL^) . These problems persist even allowing for a non- vanishing 
cosmological constant (which could resolve the age-expansion rate problem). See Figs. [LI 



and 12 



In contrast, the low-D case leads to severe tension between the SBBN prediction and the 
inferred primordial abundance of 4 He (see Figs, p], |], and [5]). This stress on SBBN can be 
relieved if the primordial helium mass fraction, derived from observations of low metallicity 
HII regions, is in error — due, perhaps, to unaccounted systematic effects — by an amount 
AF sys > 0.011 (see Fig. §). 

Alternatively, this conflict could be evidence of "new physics" |5B| (N v ^ 3; see Figs. [7| 
and §). The best fit between predictions and observations with low D is for N v = 1.9 ± 0.3. 
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One way to alter standard BBN is to change the physics of the neutrino sector. For example, 
many models predict the existence of sterile neutrinos, which interact only by mixing with 
the ordinary neutrinos. Such sterile neutrinos would not contribute significantly to the 
number of effective neutrinos (2.991 ± 0.016) inferred from the Z line-shape j5"4]| , but could 
be produced cosmologically for a wide range of masses and mixings [fjSfl . However, they only 
increase N u , exacerbating the discrepancy. 

Another possibility arises if v T has a mass in the range 10 MeV ^ M Vt < 24 MeV (the 
upper limit is the recent result from ALEPH []56|]). In this case BBN production of 4 He can 



be either increased or decreased (relative to the standard case), depending on whether v T is 
stable or unstable on nucleosynthesis time scales (~ 1 sec). An effectively stable v T (r > 10 
sec) in this mass range always increases Y relative to the standard case [[57]] (but, see p8| ) 
and would thus make for a worse fit with the data. However, if v T has a lifetime % 10 sec 
and decays into + (where is a 'majoron-like' scalar), f\ it is possible to decrease the 
predicted Y relative to the standard case (see figures 3 and 7 of Ref. [§0J)- Such an unstable 
v T contributes less than a massless neutrino species at the epoch of BBN, thereby reducing 
the yield of 4 He. For example, a v T with mass 20 — 30 MeV which decays with a lifetime of 
~ 0.1 sec reduces N u by ~ 0.5 — 1 (and Y by ~ 0.006 — 0.013, respectively), thus helping to 
resolve the apparent conflict between theory and observation. It is also possible to alter the 
yield of BBN 4 He by allowing v e to be degenerate flB~I |. If there are more v e than z/ e , Y is 



reduced relative to the standard (no degeneracy) case as the extra z/ e 's drive the neutron-to- 
proton ratio to smaller values at freeze-out. A reduction of Y of ~ 0.01 can be accomplished 
with a v e chemical potential of n e /T u ~ 0.03, corresponding to a net lepton-to-photon ratio 
of 0.005. This is to be compared to the net baryon asymmetry which is smaller by ~7 
orders of magnitude. Nevertheless, scenarios for a large lepton asymmetry are possible |p2" 



Lastly, one can relax the assumption that baryons are homogeneously distributed. However, 



'Decays with v e in the final state can directly alter the neutron-to-proton ratio and thus affect 



Yp somewhat differently |59]. 
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inhomogeneous BBN typically results in higher Yp, and therefore does not naturally resolve 
the 4 He-D conflict [|3j. 

Provided that the high-F, low-D challenge can be resolved (by AY sys > 0.011 and/or 
N u < 3), low-D is consistent with the Pop II 7 Li abundance if there has been a modest 
amount of lithium destruction/dilution in the oldest stars (see Figs. |T] and ^). The higher 
baryon density for low-D strengthens the case for dark baryons (see Fig. |9|), although that for 
non-baryonic dark matter, while still very strong, is somewhat weakened. When folded with 
the hot gas bound on the x-ray cluster baryon fraction, a "cluster baryon crisis" persists, 
arguing for Qm < 1 (see Figs. |ITiHT2l). 
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TABLES 

TABLE I. The constraints on 7/10, &>Bh 2 , N u , 4 He mass fraction (Y), and 7 Li abundance from 
the high [O] and the combined ISM and low QSO-deuterium abundances [|l3,14,2| along with those 



from solar systi 
parentheses are 


3m D and 3 He abundances Q. 
for 95% C.L. 


The errors are for 68% C.L., 


while the ranges in the 




High Dqso 


Low Dqso & Dism 


D , 3 He & BBN 


Obs. D/H (10- 


5 ) 19 ±4 


2.4 ±0.3 ±0.3, > 1.6 ±0.2 


2.6 ±0.9 


Vw 

n B h 2 

N v 
Y 

7 Li/H (10- 10 ) 


1.8 ±0.3 (1.3 - 2.7) 
0.007 ±0.001 
(0.005 - 0.010) 
2.9 ±0.3 (< 3.6) 
0.234 ± 0.002 
(0.231 - 0.239) 
1.5 ±0.6 (0.7-3.0) 


6.4±{$ (5.1 - 8.2) 
0.023 ± 0.003 
(0.019 - 0.030) 
1.9 ±0.3 (< 2.4) 
0.249 ± 0.001 
(0.246 - 0.252) 
4.7 ±0.7 (3.0 - 7.8) 


5.0 + L_| (3.5 - 7.9) 
n m 8+ - 005 

U.U1»_ 0003 
(0.013 - 0.029) 
2.1 ±0.3 (< 2.6) 

n 247 +0.003 
U-**' _o.002 

(0.243 - 0.251) 
2.9lg;g (1.4 - 7.3) 
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FIG. 1. BBN predictions (solid lines) for 4 He (Fp), D (y2p), and 7 Li (yjp) with the theoretical 
uncertainties (lex) estimated by the Monte Carlo method (dashed lines). Also shown are the 
regions constrained by the observations at 68% and 95% C.L. (shaded regions and dotted lines, 
respectively). We use the QSO measurements for y2p from Ref. [|l^] and (JUL! 
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FIG. 2. The likelihood distributions for yi implied by the low (QSO) D value (solid line) Q, 
the combined low D QSO and ISM limit (dotted line), and the combined limits from low D and 
ISM and solar system lower bounds (dashed line). 
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FIG. 4. The BBN prediction for the primordial 4 He abundance implied by each of the two QSO 
D measurements. The shaded regions and the dotted lines correspond to the 68% and 95% C.L. 
constraints. The 4 He abundance derived from HII region observations lies between the dot-dashed 
lines (68% C.L.) 
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FIG. 5. The BBN prediction for the primordial Li abundance implied by each of the two 
QSO D measurements. The shaded regions and the dotted lines correspond to the 68% and 95% 
C.L. constraints. The plateau range (see text) derived from the stars in the galactic halo (ignoring 
depletions and creations) is indicated between the dot-dashed lines. 
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FIG. 6. The two QSO deuterium constraints, combined with the BBN prediction (long-dashed 
lines, la), are shown in the U2p-y3P plane (shaded regions at 68% C.L. and dotted lines at 95% 
C.L.). The regions inside the solid, dashed, dot-dashed, and long-dashed curves are the abundances 
consistent with the solar system data for g% = 0.10, 0.15, 0.25, and 0.50, respectively ||. 
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FIG. 7. The likelihood functions for Nj, derived from the combined observations of D, 4 He, 
and 7 Li. The solid, dashed, and dotted curves are with high D, low D (QSO/ISM), and the solar 
D and 3 He ||, respectively. 
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FIG. 8. The allowed range of N u for high-D and low-D (QSO/ISM) as a function of systematic 
offsets (Ay sys ) in the 4 He abundance derived from HII region data. The shaded regions (dotted 
lines) are for 68 (95)% C.L. 
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FIG. 9. The baryon density parameter Ob versus the Hubble parameter Hq. The two bands 
correspond to the 68% (shaded) and 95% (dotted) C.L. ranges for r/io inferred from SBBN for high 
Dqso and low Dqso/Dism (see Table I). Also shown are the estimates [|36| of the contributions to 
Qb from luminous baryons in galaxies (dashed curve) and a dynamical estimate |37]] of the lower 
bound to the total mass density parameter (solid line). 



26 



1.0 




FIG. 10. BBN and x-ray cluster constraints on the total (matter) mass density parameter 
(Om) versus the Hubble parameter (Hq) for the two choices of primordial D. The shaded bands 
(dotted curves) are the 68% (95%) C.L. allowed regions (upper limits; see text). 
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FIG. 11. The x-ray cluster constraints on the total mass density parameter (Qm) versus the 
Hubble parameter (Hq) for high and low D (shaded bands, 68% C.L.) along with the constraints 
from the shape parameter T = £ly[h and the lower bound to Qm from dynamics (Odyn)- Also 
shown are the Om vs. Hq relations for two choices of the present age of the universe. 
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FIG. 12. As in Fig. 11, but for a zero curvature (k = 0) model with A ^ (Om + = 1) 
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